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I Abstract. This is a continuation of the work initiated in [18j on so-called driven 

1^ ' cofactor systems, which are partially decoupling second-order differential equations 

^ • of a special kind. The main purpose in [T^ was to obtain an intrinsic, geometrical 

G . characterization of such systems, and to explain the basic underlying concepts in a 

i brief note. In the present paper we address the more intricate part of the theory. It 

I involves in the first place understanding all details of an algorithmic construction of 

' quadratic integrals and their involutivity. It secondly requires explaining the subtle 

■ way in which suitably constructed canonical transformations reduce the Hamilton- 

. Jacobi problem of the (a priori time-dependent) driven part of the system into that 

! of an equivalent autonomous system of Stackel type. 

o ; 
^ , 

^ ; 1 Introduction 

^ I To the best of our knowledge, the idea of a cofactor system stems from a paper by Rauch- 
■ — ' Wojciechowski et al. on certain Newtonian systems in the Euchdean plane [16], but the 
term 'quadratic integral of cofactor type' came in use with Lundmark's generalization to 
systems of arbitrary dimension, which was published with some delay in [5]. The term 
refers to the fact that the matrix of the quadratic part of the first integral comes from the 
cofactor tensor of a tensor which has special properties with respect to the (Euclidean) 
metric of the kinetic energy of the system (a tensor which was termed 'inertia tensor' 
in Benenti's work [1] on Stackel systems). The other point to be emphasized about the 
'quadratic integral of cofactor type' is that the zeroth-order terms in the integral do not 
come from a potential energy function, in other words, the Newtonian systems under 
consideration have force terms which are of nonconservative type, albeit of a very special 
nature, determined by a scalar function also and by the cofactor tensor. It was recognized 
in [5] that the work of Lundmark could easily be generalized to kinetic energy terms 
associated to an arbitrary Riemannian metric and in fact fits perfectly within the theory 
of 'special conformal Killing tensors', as studied in [B]. This gave rise to an intrinsic 
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characterization of general cofactor systems, in the context of which it was shown that 
the corresponding quadratic integral in fact is the Hamiltonian for a quasi-Hamiltonian 
representation of the system with respect to a non-standard Poisson structure, coming 
from the special conformal Killing tensor (scKt for short). It must be said that scKts for 
general metrics had already appeared also, with slightly different assumptions and under 
different names, in the work of Benenti (see the comprehensive review paper |3]). 

As our title reflects, there is an extra aspect about the cofactor systems we want to 
study. So called 'driven cofactor systems' again were introduced by Lundmark and Rauch- 
Wojciechowski [10], still in the context of mechanical systems with Euclidean kinetic 
energy metric and hence having no terms quadratic in the velocities in the second-order 
equations of motion. Briefly, the systems discussed in [10] are of the form 

= Q'iy^), i = l,...,m 
F = Q°(?/\x^) a = l,...,n. 

To begin with, they exhibit a given partial decoupling whereby the ^/-equations are referred 
to as the driving system, and the x-equations as the driven part. In addition, it is 
assumed that the overall system is of cofactor type and that, more restrictively, the force 
terms come from a potential, parametrically depending on the driving coordinates. 
These are rather strong conditions indeed, but they were shown to lead to quite striking 
conclusions in [10]. First of all the driving system turns out to be of cofactor type in 
its own right. Secondly, the driven system, when regarded as a time-dependent system 
along solutions y(t) of the driving system, has n (time-dependent) quadratic integrals. 
Most astonishingly, however, the authors managed to show that (under some technical 
assumptions) there exists a time-dependent standard canonical transformation, which 
has the effect of shifting the time-dependence to an overall factor, so that an autonomous 
Hamiltonian can be identified which turns out to be of Stackel type. 

Clearly, such results lead to a double challenge. The first question is whether they can be 
extended to general cofactor systems and can be understood in more intrinsic geometric 
terms this way. Secondly, if there exists an intrinsic scheme behind these observations, 
coordinates with respect to which the system partially decouples should not be regarded 
as part of the data; their existence and a constructive procedure to find them should follow 
from testing coordinate free conditions. This second aspect fortunately has sufficiently 
been explored in the literature. Second-order equations with such a decoupling property 
were called submersive in [8], where both local and global criteria were developed for 
their characterization. The local conditions for existence of appropriate coordinates and 
for their construction were turned into a more compact and transparent form by using the 
geometric calculus along the tangent bundle projection developed in [13] (see also [I2|)- 

In a brief communication [TB|, one of us has presented the main ingredients towards the 
resolution of the above double challenge. First of all, [IB] provides a coordinate free 
definition of a driven cofactor system. It further shows that the driving system carries a 
cofactor structure in its own right and that the overall system has a second, in some sense 
degenerate cofactor representation, from which one can develop a scheme to construct n+1 
quadratic integrals. One of these is the integral corresponding to the cofactor nature of 
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the driving system; the other n (the dimension of the driven system) are integrals of the 
driven system along solutions of the driving one. What was not sufficiently understood 
at that time was the delicate issue of the nature of a canonical transformation (if any) 
which would have the effect of eliminating the time-dependence form the driven system 
in a way which leads to the identification of a Stackel system, as in the Euclidean case. 
The main purpose of the present paper is precisely to complete this part of the story. 
To understand precisely what is happening in the process, we need much more explicit 
information about the recursive procedure which leads to the n + 1 first integrals referred 
to before, and also about the precise structure of these integrals. The details of such 
computations will take a substantial part of our discussion, but in doing so, the present 
paper will also complement a number of the results already mentioned in [18] , where the 
focus was more on existence issues. 

The next section contains the basic definitions needed to describe driven cofactor systems 
and a summary of results presented in [18]. In Section 3, we develop the algorithm which 
leads to the identification of the n + 1 first integrals referred to before. A key issue for 
understanding the nature of the driven system is the identification of a scKt for its proper 
metric. It is shown in Section 4 that this scKt does not, however, turn the driven system 
into a cofactor system in its own right; instead it does that for a modified driven system 
which turns out to play a role further on. In Section 5, we start by identifying Darboux 
coordinates for the symplectic form associated to the scKt of the complete system, they 
are obtained by suitably modifying the momenta. But we gradually develop arguments 
then to come to an even better selection of modified momenta, which takes the specific 
decoupling properties of our system into account and are shown to be related to a time- 
dependent (standard) canonical transformation for the driven part of the system. In 
Section 6, we prove that the application of this canonical transformation, followed by one 
which comes from using eigenfunctions as new coordinates, produces the rather miraculous 
effect of reducing the driven system essentially to an autonomous Stackel type system. 
The proofs in Section 6 are partly based on simple, indirect arguments, but they are 
supported also by explicit computations about the structure of all first integrals, which 
are presented in an Appendix. A couple of illustrative examples are presented in Section 7. 



2 Preliminaries 

The mechanical systems we are talking about in this paper belong to the class of noncon- 
servative Lagrangian systems, governed by equations of the form 

"^^^^ (1) 



dt \dq°'J dq' 

Here T = ^gap(q)v"'v'^ is a kinetic energy function on the tangent bundle TM of a 
Riemannian manifold M with metric g; the nonconservative forces Qa are assumed to 
depend on the position variables only and thus can be viewed as components of a 1-form 
Ai = Qa{,q)dq°' on M. We recall from [6] that a special conformal Killing tensor J (scKt) 
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is a type (1, 1) tensor field on M, which is symmetric with respect to g and satisfies 
(lowering an index in the usual way) 

<^a/3|7 = ^{c^aQiS"/ + ct/sQa-y), which further implies that a = diti J). (2) 

The Nijenhuis torsion A/j vanishes, implying that dj'^ = 0. Recall for completeness (see 
[7j) that dj is the derivation of degree 1 which (anti) commutes with the exterior derivative 
and whose action on functions / G C°°(M) is given by djf = J{df). The Nijenhuis torsion 
A/j is the vector-valued 2-form with components 

ja,dJl_djf\_j,dJl^j^dJ^ 



We also have that Dj^ = 0, where Dj is the 'gauged differential operator', defined by 

Djp = djp + d{tT J) Ap = (det J)~Mj((det J)p). (3) 
The equality of both expressions in the defining relation of Dj follows from the fact that 

(ij(det J) = (det J)d(tr J), (4) 
for any tensor with vanishing Nijenhuis torsion. 

The concept of a cofactor system on a general Riemannian manifold was introduced in [5] 
in the following way. 

Definition 1. A cofactor system is a triple {g,p, J) on a Riemannian manifold where g 
is the metric, p is a 1-form and J is a nonsingular special conformal Killing tensor such 
that Djp = 0. 



Combining this concept with the geometric notion of submersiveness, we came in [18] to 
the following generalization and coordinate free formulation of the kind of driven cofactor 
systems introduced in [TO] . 



Definition 2. A driven cofactor system is a cofactor system {g,p,J), for which there 
exists a distribution K along the projection r : TM — >■ M , with the properties 

^{K) CK, VK C K, D'zK CK VZ G X{t), (5) 

dp{K, K) = 0, DXA'^, K) ^ 0. (6) 

To understand the meaning of these conditions, one needs to know about the following in- 
trinsic geometrical concepts associated to general second-order differential equation fields 
(soDEs), say 

Such a SODE comes with a connection, determined by the following horizontal lift con- 
struction from X[M) (the module of vector fields on M) to X(TM) 
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where 



H^ = ^ r^^^, with n 



are the connection coefficients. Naturally, we also have a vertical lift 
X = = X° Va, where 



It is clear that both operations still make sense if we allow the components of X to be 
functions X"(g, v), meaning that we extend the domain of the horizontal and vertical lift 
to A:'(t), the C°°(TM)-module of vector fields along r : TM — )■ M. These lifts further 
give rise to corresponding horizontal and vertical covariant derivative operators and 
Dx, determined by the following action on functions F G C'^{TM) and basic vector fields 
(and then further extended by duality): 

D-F = X"K(F), D-A = o 

The in ([6]) is a covariant differential, defined for any tensor field T along r by 
D"T{X, ¥,...) = DxT{Y, . . .). Furthermore, the decomposition of CrX'^ into its horizon- 
tal and vertical part identifies the important concepts of dynamical covariant derivative 
V, a self-dual degree derivation on tensor fields along r, and Jacobi endomorphism, a 
(1,1) tensor $ along r: 

£rX^ = (VX)" + $(X)^. 
For practical purposes, it suffices to know that: 



For the broader picture of derivations of forms along r, where the above mentioned con- 
cepts play a distinctive role, one can consult [HI [T4] . 

As was discussed in Proposition 1 of [18] and before that already in [12], the existence of 
a distribution K along r having the properties ([S]) implies that K is actually spanned by 
a distribution on M which is integrable, and introducing adapted coordinates (?/*, x") has 
the effect that the given SODE partially decouples into equations of the form 

f = r{y.y). z = i,...,m, (7) 

-a _ f"-(^x,y,x,y), a = l,...,n (here n + m = X = dim M). (8) 

We will keep referring to the decoupled ^/-system as the driving equations and to the 
remaining x-equations as the driven system. But there is more to it in this context, which 
brings us to the final requirements (E]) in definition 2. of course is the orthogonal 
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complement of K with respect to the Riemannian metric g. An important property of g 
is that 

Vg = 0, (9) 

which follows from the fact that g is the Hessian of the Lagrangian for the conservative 
part of the system and the knowledge that the nonconservative forces determined by 
have no effect on the connection coefficients. It was shown in [TB] that as a result of 
([9]) K-^ inherits the properties VK-^ C K-^ and D^/^-*- C K-^ from K and that this is 
enough to conclude that the two complementary distributions are in fact simultaneously 
integrable. Hence we can choose x°' and coordinates which are adapted to K and K-^ 
at the same time, i.e. in such a way that 

It also follows that the kinetic energy part in the equations of motion ([1]) decouples com- 
pletely. In other words, if we put gi := g\j^i_ and g2 := g\j^, then in adapted coordinates: 

9i = 9ij{y)dy' ® dy\ g2 = gab{x)dx'' (g) c^a;^ (10) 

while gia = gai = 0. Similarly, for the corresponding connection coefficients, we have 

r}, = r},(y), Tl = TUx) 

and all other connection coefficients are zero. If we insist on keeping some partial coupling 
between a driving and a driven part of the dynamics therefore, this can only come from 
the nonconservative forces in /i. The second of the conditions ([6]) exactly guarantees such 
a coupling. The condition dfi{K, K) = on the other hand can easily be seen to model 
the additional assumption that, in adapted coordinates, the driven part has force terms 
Qa which are derivable from a potential energy function (with parametric dependence on 
the driving coordinates 



3 The cofactor pair scheme on M and n + 1 quadratic 
integrals 

It is appropriate to introduce complementary projection operators 

Pi : X{M) -> K^, P2 : X{M) K. 

We thus have Pi + P2 = In (the identity tensor on the A^- dimensional manifold M) , 
Pi o P2 = P2 o Pi = 0, P^ = Pi, and we occasionally put Pi\k^ = Im, P2\k = In- As in 
[18] , we then look at the scKt tensor J (for its action on vector fields) as the sum of the 
following four parts 

J, = PiOjoP,, 1=1,2, J12 = Pi O J O P2, J2l=P20j0P,, (11) 
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and we shall use a similar notation also for other type (1,1) tensor fields of interest further 
on. 



Important notational convention: In principle, such tensor fields act on the whole 
module of vector fields on M; we shall use the same notation, however, when we consider 
their restriction to the appropriate submodule K or K-^ where they are not zero. 

Recall that the scKt conditions (E]), when expressed in terms of the coefficients of the 
original type (1, 1) tensor, take the form 



QJ<^ 

Jp\, ■■=^- J^n-y + JF''^, = |(«/3^7 + <^^9'"'9^,)- (12) 



Taking the decoupling properties of g into account, it follows that in adapted coordinates: 
dJi/dx'' = dJ^/dy' = and 



81°^ dp 

J^\^ = ^- JFlk = ¥y>9'"'9^.. ■% = ^a- JT,a = h^9''9a,. (13) 



Hence, the different blocks of J have the following type of restricted dependence on the 
adapted coordinates: 



while 



with 



d ■ ■ d 

J21 = J^iy,x)— ® dy\ J12 = Jl{y,x)—^ ® dx\ 



^ . and — = — ^. (14) 

Furthermore, with /ij := Pi(/i) and assuming that Ji is nonsingular, we know from [18] 
that (5'1,/ii, Ji) provides a cofactor representation of the driving system. This implies in 
particular (see further for more details) that this system has a quadratic integral 

E' = \Al{y)yY + W\y), (15) 

where = cof Ji. For clarity: the cofactor tensor A of a type (1, 1) tensor J (notation 
A = cof J) is defined by the relation J A = AJ = (det J)I. The function is determined 
by the relation A^fii = —dW^, which is locally equivalent to the condition Dj^fii = in 
the definition of a cofactor system. Incidentally, the minus sign in the expression for A^fii 
is in correspondence with the in the integral (and [18] contains a sign error in 

this sense). 

We finally recall from [18] that P2, being a degenerate type (1, 1) tensor field on M, 
formally satisfies the requirements for a scKt with respect to g, and additionally has the 
property Dp^fi = 0, meaning that we have a second (degenerate) cofactor representation 
for the full nonconservative system F, the implications of which we will investigate now. 
The starting point is that J + aP2 also satisfies the scKt condition for any constant a (and 
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is nonsingular for sufficiently small values of a). Let A{a) denote the cofactor tensor of 
J + aP2, so that 

(J + aP2)A{a) = A{a){J + aPa) = det(J + aP2)lN- (16) 

Since P2 = In in adapted coordinates, it is clear that A{a) and det( J+aP2) are polynomials 
in a of degree n. We represent them as follows, 

n+l n+1 

A{a) = ^ A(i)a*"\ det( J + aPa) = ^ A(i)a*"\ 
and identify the coefficients of equal powers of a in the identities (fT6|) . We get 

P2A(„+i) = A(„+i)P2 = 0, (17) 

JA(i+i) + P2v4(i) = A(,+i)J + A(i)P2 = A(,+i)/7v, il<i< n), (18) 
JA(i) = A(i)J = A(i)/;v (19) 

Information about the block structure of the different should follow by left and right 
actions of the projectors P^ on these relations. It immediately follows from (1171) for 
example that 

A(n+l)21 = ^(n+l)12 = ^(n.+l)2 = 0, (20) 

while (ITSl) with i = n subsequently implies that Jiy4(„+i)i = A(^n+i)iJi = A(„+i)Pi, or in 
the restriction to K^: 

Taking into account that terms of degree n in A{a) can only be produced by minors of 
the Ji elements of J, this implies that 

A(n+i)i = A^ = cof Ji and A(„+i) = det Ji. (21) 

The same equation flTS]) for i = n then further yields information about parts of it is 
indeed not hard to show by appropriate actions of the projectors that (in the restriction 
to K or K^) 

A{n)2 = (det Ji)J„, A(„)2i = -J2i(cof Ji), A(„)i2 = -(cof Jl) Jl2. (22) 

For the remaining part of we have to move to the next line in the hierarchy (i = n — 1 
in f|T8l) ): a right and left action of Pi leads to JiA(^n)i + ^i2^(r4)2i = A(„)Pi, from which 
it readily follows that 

= (det Ji)J{'Ji2J2iJ{' + A(„) (23) 

We will come back in more detail to the continuation of this recursive scheme when we 
are in a position to gather information about the functions A(j). But note that the final 
relation ( fT9l) expresses that = cof J and A(i) = det J. 

The important feature about having a cofactor pair scheme is that it gives rise to a 
family of quadratic integrals which are in involution with respect to a double Poisson 
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structure. We sketch how this will work here by following the procedure explained in [5]. 
Let g : TM T*M denote the Legendre type diffeomorphism coming from the given 
metric g, which in (general) coordinates reads pa = gapv^ , and put V = ^*r, where T is 
the SODE defined by the cofactor system {g, fi, J). Then, according to theorem 2 in the 
requirements about a cofactor system are equivalent to T representing a quasi-Hamiltonian 
system with a quadratic Hamiltonian with respect to the Poisson tensor which can be 
constructed from the tensor J with vanishing Nijenhuis torsion. More precisely, if J 
denotes the so-called complete lift of J to T*M and Pq is the canonical Poisson map on 
T*M, i.e. the canonical Poisson tensor, interpreted as a map Pq • X*(T*M) -> X{T*M), 
we put Pj = J o Pq and the point about T is that 

(det J)f = Pj{dH), with H = lA'^^p^pp + W, (24) 

where A is the cofactor tensor of J and W is the function, determined (up to a constant) by 
the property A^ = —dW. Having a double cofactor representation with special conformal 
Kilhng tensors such as J and P2 in our present situation, will entail that V satisfies a 
relation of the form 

(det(J + aP2))V = Pj+aP2dH{a), with H{a) = lA'^^{a)p^p(s + W{a). (25) 

Since Pj+aP2 = Pj + (^Pp2^ one can see that H{a) will be a polynomial of degree at most 
n in the parameter a, say of the form 

n+l 

H{a) = Y,H^i)a'~\ with = + (26) 

i=l 

where comes from the A{i) tensor considered before, with an index raised by the given 
metric g, and W(i) is determined by the relation = —dW{i). Naturally H{a) and 

therefore all its coefficients H(i) will be first integrals of the system f . 

Before proceeding, observe that H(n+i) is essentially the function which is a first 
integral of the driving system (more precisely -ff(n+i) = g*E^). Indeed, it follows from 
the block structure of A^n+i) determined before that in adapted coordinates, with pi = 

= l{A'f{y)p,p, + W\y). (27) 

It is worth mentioning also that for a general J, the coordinate expression of Pj is given 
by 

' ^ dpp dq- ^ J gp^ gp^ ■ 



Naturally, in adapted coordinates, 

^ d d 



P2 



dpa dx°- 

where the momentum variables Pa are defined by Pa = gab{x)x^- 



A 71-7, (29) 
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Taking the various polynomial representations into account, property (l25l) becomes 

(Er=YAw«^"')f = {Pj + aPpj(Er=i'^^»«^"')- 

Identifying coefficients of equal powers of a requires ffist of all that we should have 
Pp^{dH(^n+i)) = 0, and this is clearly verified in view of the preceding observations. We 
further must have that 

A(,)f = Pj{dH(i)) + Pp,{dH^,_,)), 1< ^ < n + 1, (30) 

and finally for i = 1 that A(i)r = Pj{dH(^i)), but this is merely a conffimation of the 
quasi-Hamiltonian structure coming from J, since A(i) = det J and //(i) = H. There is 
no reason to expect that the n ffist integrals up to of T would depend on the 
coordinates {x°',pa) only. Nevertheless, they will be ffist integrals of the driven system 
along solutions {y'' (t) , Pi(t)) of the driving system. As for the question of involutiveness, 
if we adopt the notational convention that Pj{df) = {/, ■ }j, it follows from (1301) that 

{H^r),H^i)}j + {H^i^i), H^i)}p, = A(i)f(i7(o) = 0, l<z<n + l,l</<n + l, (31) 

and from the two other observations about and -ff(i) that 

= {H^n+i),H^i)}p, = V/. (32) 

Using ( l32i) . it further follows from ( 13T]) with / = 1 and I = n + 1 respectively, that also 

{iJ(i),if(i)}p2 = {H^i),H^n+i)}j = 0, 

and then a simple recursive argument in ( l3T]) finally implies that all are in involution 
with respect to both the J-bracket and the P2-bracket. It is worth observing (see ( 129|) ) that 
in adapted coordinates the P2-bracket formally looks like the standard Poisson bracket in 
the {x"',pa) coordinates so that we have, along solutions of the driving system, n integrals 
for the driven system which are in involution in the standard sense. 

4 An induced special conformal Killing tensor for the 
metric of the driven system 

Part of our basic assumptions so far is that both J and Ji are nonsingular. We shall see 
now that this naturally leads to the introduction of another nonsingular type (1,1) tensor 
J2, which is a kind of deformation of J2. When talking about nonsingularity here, don't 
forget the notational convention specified before! To say that Ji is nonsingular of course 
only makes sense when we mean Ji\k-^', likewise, nonsingularity of J2 will refer to J2\k- 

For later reference we look at the adjoint action of J on 1-forms for a moment. In doing 
so (as we already tacitly did for the cofactor tensor A) we use the same notation again, 
i.e. do not write J* as is sometimes customary. One has to keep in mind, however, 
that when compositions are involved (as in the definition of J12 and J21), the order of 
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such compositions has to be reversed. Let a be an arbitrary 1-form on M and put 
(3 = Ja = J{Pia + P2a). To solve such a relation for a in terms of (3, it is natural to 
actually compute Pia and P2«- We have 

Pi/3 = Ji(Pia) + J2i(P2«), 
P2/3 = MPia) + J2{P2a). 

Since Ji is nonsingular, it follows from the first relation that 

Pitt = J^\Pil3) - Jr'</2i(P2«), (33) 

and substitution of this result in the second relation leads to 

J2(P2«) = (P2 - Jl2Jl'Pl)/3, (34) 

where J2 is defined (for its action on 1-forms) as 

J2 = J2 — J12J1 ^^21- (35) 

Obviously J2 vanishes on i^-*-, but will be nonsingular on K so that P2a now can be 
obtained from flM|) and substitution in fl33|) subsequently gives us Pia. In fact, what we 
are looking at here is the following typical factorization of a matrix with a block structure, 
this time written as representing the action of a type (1, 1) tensor on vector fields: 

f Ji Ji2\^fJi \ / 1 JfVia^ 
\ J21 ^2 ) \ 'hi 1 / V J2 — J21J1 J12 

It follows that det J = (det Ji)(det J2) so that nonsingularity of J and Ji implies the same 
for J2. In adapted coordinates, J2 has components 

j2:i = j^-jnJi'pl- (37) 

Proposition 1. J2 is a (parameter dependent) special conformal Killing tensor for g2, 
and its cofactor tensor is (det Ji)~^A2 with A = cof J. 

Proof. The proof is a straightforward computation, for which it will be suitable to work 
in the adapted {y\x°') coordinates. Let us lower an index to apply the scKt condition in 
the covariant form ([2]). Keeping in mind that Qai = g-ia = 0, we have 

hch ■= Qcahh — Jch ~ Jci{Jl "^Y^ Jjb- 

Since J is symmetric, the same is true for J2. Note that J2 depends on both sets of 
coordinates, but the are regarded as external parameters for our present considerations. 
Since Ji depends on the only and the same of course holds for its inverse (or its cofactor 
tensor A(n^i)i), we get in the first place that 

J2cb\a — Jcb\a Jci\a{Jl Y'^'^jb Jcii'Jl Jjb\a- 
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It follows from (|2]) that 

Jcb\a = lioicOba + Oibgca), and Jci\a = ^OiiQca- (38) 

We then easily obtain that 

J2cb\a = \{acgba + "bfi^ca), with = - JciiJl^Y^Oij, (39) 

which shows that J2 is a scKt for g2. 

If we denote cof J as in (!24|) by A (and we have seen that it is also ^(1)), we know that 
(det J)g'''' = J^A^"" = J^A'" + J^A''" and = (det J)g^'' = J^A''^ + JlA^^. It then follows 
that 

hlA'^detJ^r^ = {deiJrr^{JiA''-Jt{J^%JiA'') 

= (det Ji)-i ((det J)(7"^ - J^A'^ + Jf ( Jf i)} J^A'^'^) 
= (det Ji)-^(det J)^'^'^ 
= (detJs)^?"^ 

which proves the last statement of the proposition. □ 

Note that J2, perhaps rather unexpectedly, does not give rise to a cof actor system repre- 
sentation of the driven system in a strict sense. In other words, it does not seem to be true 
that the forces 112 of the driven system have the property Dj^^2 = 0, or equivalently that 
(cof J2)/^2 is closed. At this moment, the closest we can get to such a property is that an 
associated driven cofactor system {g2,p'2, J2) exists, with modified nonconservative forces 

Proposition 2. A cofactor system, parametrically depending on the coordinates of the 
driving system, is determined by the triple {g2,P'2, J2), where g2 and J2 are as before, and 

Ji2 = P2{dWi^n)). (40) 

W(n) being the function encountered in the recursive scheme following from ( f^) . 

Proof. We know that the given nonconservative forces yU. satisfy the relation A^ = —dW, 
or equivalently (det J)fi = —J{dW). Consider the scheme of the beginning of this section 
which led to the introduction of J2. Letting —dW play the role of a and (det J)/i the role 
of /3, we put Pi{dW) = diW for convenience. The relation becomes 

- J2{d2W) = (det J)/i2 + (det J2) MdiW^) , (41) 

where for the second term on the right we have taken into account that fii satisfies 
(det Ji)fii = —Ji{diW^) in view of the cofactor representation of the driving system, and 
that det J = (det Ji)(det J2). Secondly, projecting the relation = —dW(n) under P2, 

we get in the first place that 

^(n)2/^2 + ^(n)12/^l = " C^2W^(n)- 
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Using the information gathered about in fl22|) and the cofactor system property of 
the driving system which was just recalled, this immediately leads to 

(det Ji)/i2 + MdiW^) = -d2W^n). (42) 

Substituting this result in ( HTj) . we obtain the relation 

{detJ2)d2W^n) = J2{d2W), (43) 

which implies, with ^2 = d2W(^n), that 

(Cof J2)fi2 = d2W. (44) 

Together with the knowledge that J2 is a scKt with respect to g2-, this expresses that the 
triple {g2,p'2, J2) determines a cofactor system. □ 

It is worth emphasizing again that the cofactor representation of this modified driven 
system must be seen also as a statement about second-order differential equations for the 
which parametrically depend on the |/-coordinates. This is clear, for example, from 
the fact that the exterior derivative in the picture is ^2- 

As explained in the introduction, one of our main objectives is to explore the remarkable 
situation that the driven system, although being essentially time-dependent along solu- 
tions of the driving system, does give rise in the end to a Stackel-type Hamilton- Jacobi 
separability anyway. This will require a supplementary assumption about existence of 
independent eigenfunctions of J. We want to show at this point that it is again the tensor 
J2 which is relevant for this purpose. 

Consider the (degenerate kind of) eigenvalue equation det(J — AP2) = 0, which is a 
polynomial equation of degree n for A. Our basic assumption now is that this equation 
has n functionally independent solutions in the sense that the 1-forms P2{du"') are 
linearly independent. If we think of the m° as expressed in terms of the adapted {y^,x^) 
coordinates, this amounts to saying that the Jacobian {du°'/dx^) is nonsingular. From 
the identity det( J — u°-{y, x)P2) = for each fixed u"(y, x), it follows that 



= rf(det(J-u"(?/,x)P2)) 

(9(det(J- AP2)) 



d{dei{J - \P2))\x=u-{y,x) 



du". 

\=u°-(y,x) 



dX 

But since J — XP2 has vanishing Nijenhuis torsion, we know from (jlj) that 

(J - AP2) d (det(J - AP2)) = det(J - AP2) dtr(J - AP2) 
for all A. It then follows by acting with (J — u"-{y, x)P2) on the preceding identity that 
a(det(J- AP2)) 



(J-M"(y,x)P2)rfM" = 0, 

\=u'^(y,x) 



d\ 

and thus, since all eigenfunctions are assumed to be simple, that 

{J -u''{y,x)P2)du'' = 0. (45) 
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Proposition 3. Assume that the equation det(J — XP2) = has n functionally indepen- 
dent eigenf unctions u"" . Then du°' is an eigenform of J (in the sense of equation [J^)) 
corresponding to the eigenvalue u"". Moreover, the u"" are also eigenf unctions of J2, with 
P2{du°') as corresponding eigenform. 

Proof. It remains to prove the statement about J2. For that, it suffices to go back once 
more to the analysis about /3 = Ja at the beginning of this section, with du"" in the role 
of a and u°'P2{du°') (no sum!) in the role of (3. The relation (134|) then says that 

J2{P2{du'')) = U»P2(rf^"), 

which is precisely what we need. □ 



5 A symplectic view-point and Darboux coordinates 

Since J is assumed to be nonsingular, the Poisson tensor associated to Pj actually comes 
from a symplectic form which we call uj. The sign convention which we adopt here is 
that for any function F 

X = Pj{dF) ^ ixooj = -dF 

One easily verifies that, referring to the general coordinate expression (l28l) of Pj, ujj is 
given by 

W~d^) ^ ^^^^ 

So far, this correspondence is valid for any nonsingular type (1,1) tensor J on M. The 
ffist term in ujj strongly suggests introducing new momentum variables pa by 

Pa = J^{q)Pp- 

Lemma 1. The coordinate change {q,p) ^ {q,p) determines a Darboux chart for uj if 
and only if Mj = 0. 

Proof. From = J~^pPa, it is easy to compute that 

dp, A dq^ = J-';dpa A dq^ - Ip, {^J^'l " ^q" A dq" . 

Subtracting this from (H6i) . one easily obtains that the coefficient of the resulting 2- form 
(after multiplication by two J-factors) will be zero if and only if Mj = 0. □ 

As explained in Section 3, the fact that the SODE F on TM satisfies the requirements of a 
cofactor system is equivalent to saying that its image F under the Legendre map g has a 
quasi-Hamiltonian representation ( 124|) with respect to Pj. This in turn translates within 
the symplectic view-point to (det J)ipuj = —dH. The preceding lemma then says that 
uj will take the form of the standard symplectic form on T*M when expressed in the 
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variables {q,p). However, it is not clear that we can take advantage of such a coordinate 
change because it does not take account of the special feature of partial decoupling which 
our system exhibits. We shall show that there is a better choice of new momenta, which is 
inspired by the above transition to Darboux coordinates but does take the extra features 
of a driven cofactor system into account. Note in passing that there exists a different 
technique for obtaining a Hamiltonian representation out of a quasi-Hamiltonian one: 
roughly it consists of absorbing the overall factor by a change of timescale, which has the 
disadvantage, however, that the definition of the new time makes sense only along (the 
unknown) solutions of the system. This technique is well documented in [3], for example, 
and was extensively used in the context of cofactor systems in [H]. We believe that the 
line of approach we adopt here offers more insight in understanding the delicate aspects 
of the driven nature of our cofactor system. 

Going back to the SODE F on TM, we can first pass to the coordinates (y*, x") adapted 
to the complementary distributions K-^ and K, before passing to the quasi-Hamiltonian 
representation T. The projectors Pi and P2 have corresponding actions on TM through 
their complete lifts; they give rise to a partial splitting of F in the form F = Fi + F2 say, as 
exhibited in the equations ([71 [8]). Likewise, the complete lifts Pi, P2 of the projectors to 
T*M produce a partial decoupling F = ri-|-r2, which in adapted coordinates is simply the 
effect of transforming ([71 [8]) to equivalent first-order equations by passing to the momenta 
Pi = gij{y)y^ and pa = gab{x)x^. It is worth illustrating this in more detail as follows. The 
SODE F associated to the equations of motion of the form ([1]), after applying the overall 
Legendre map g, transforms to 

f = *"%|j + r;.J"W.|; + Q»|;. (47) 

In adapted {y, x)-coordinates, in view of the way the components of g and the connection 
coefficients decouple (see (fTOl) and its consequences), this expression becomes 

ax" Opb Opb 

The first line reflects the fact that the driving system has its own cofactor representation, 
i.e. satisfies 

(det Ji)fi = Pj,(rf/7(„+i)), (49) 

with -ff(n+i) as in (1271) . Concerning the second line, we should take into account the 
extra assumption that the driven system has a standard Hamiltonian representation: as 
indicated before, the condition dfi{K, K) = in definition 2 expresses that the force 
terms Qa of the driven system are derivable from a potential energy function V{x, y) say, 
depending parametrically on the driving coordinates It is then clear that the second 
line simply expresses that 

T2 = Pp,{dh), with h=lg''\x)paPb + V{x,y), (50) 
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keeping in mind that Pp^ in adapted coordinates merely is the standard Poisson structure 
in the variables (x",pa)- With this splitting of F in mind, it looks more appropriate not 
to spoil the decoupled feature of the driving system by introducing Darboux coordinates 
for the overall symplectic structure ujj. Instead, we can put Pi = Ji^Pj, which will have 
the effect of introducing Darboux coordinates for wj^. As for the driven part let us 
first investigate in detail what the introduction of the momenta p would do. 

Formally we can regard the transformation formulas Pa = Ja{Q)Pi3 &s representing a 
relation between 1-forms on M, of the type /3 = Ja discussed at the beginning of Section 4. 
It then follows from the considerations leading to ( !34l) that 

^2aPfe =Pa- JUl^lPj- (51) 

This suggests that the more relevant momentum variables for the driven system actually 
are pa ■= J-iaPb- The conclusion from this preliminary analysis is that we shall consider 
the following linear change of momenta 

P^ = J^Pj, (52) 
Pa = Pa + JaPi- (53) 

It turns out that the transition from pa to pa, viewed as time-dependent transformation 
along solutions of the driving system, actually represents a time-dependent canonical 
transformation for the driven system in the standard sense and hence is ideally suited 
to preserve the special assumption on that system. Indeed, in view of the second of the 
properties (|T^ . we know that the components J* of J can be written as J* = dip^/dx"" 
for some functions tlj^{x,y). Defining F{x,p,t) by 

F{x,p,t) = x^pa + ^ix^t), with ^{x,t) =tlj'{x,y{t))pi{t), (54) 

we create a generating function of mixed type (depending on the old position variables x 
and new momenta p) for a standard canonical transformation {x°',pa) f-)- {x°',Pa), which 
does not change the coordinates, transforms the momenta according to (jSSD, but must be 
viewed as time-dependent along solutions of the driving system. The Hamiltonian of the 
transformed system then is given by 

OF 

h{x,p,t):=h + —. (55) 

We shall see in the next section that this canonical transformation is one of two steps 
which are required to relate the original Hamiltonian h of the driven system to the first 
integral of T and that is the key to understanding the subtle way in which the 
driven system in the end corresponds to an autonomous Hamiltonian system which is 
separable in the Hamilton-Jacobi sense. 
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6 Separability of the Hamilton-Jacobi equation for 
the driven system 



Let us start by computing the function expressed in the variables {y\ x°',Pi,Pa). We 
have 

From ( l22ll . raising an index, we learn that 

Al^) = (det J,)g^^ A^^) = -A'lf-. 
Making the substitutions (152 | l53|) it then readily follows (remember that A^ is cof Ji) that 

H^n) = |(det Ji) g^^PaPb - |(det Jl)J^'4p^PJ + \Af^-^J^,^^piPj + 
Using (1251) now, we arrive at the following result: 

= |(det Ji) + \^in) f'PiPj + IV(n). (56) 

So, introducing the new variables has the interesting effect of eliminating the terms in 
mixed momenta in H(n). Obviously, however, the effect on h will be the opposite. We get 

h = Ig^'PaPb + r'PaP^ + J^P^pJ + V. (57) 

But for the interpretation as time-dependent canonical transformation, we need to look 
at the function h and will show now that this function is more closely related to -ff(„,). 

Lemma 2. Under the canonical transformation with generating function (5^ , the trans- 
formed Hamiltonian h of the driven system takes the following form, to within an additive 
function of time, 

h = (det JlY^H^n) + r'paPi- (58) 



Proof. We need to add to the expression fl57l) for h the term 

dF .f, 



dt dt dy^ 

computed along solutions of the driving equations. Since the pi were introduced to provide 
Darboux coordinates, the Poisson map Pj^ in (149|) takes the form of the standard Poisson 
map so that the equations of the driving system become 

detJi 2/ = — , detJi Pfc = 

opk oy" 

whereby the function H{.n+i) from ( l27j) . when expressed in the new momenta, reads 

H^n+i) = Kdet Ji)Jl' p^p^ + W\ (59) 
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It is now fairly straightforward to compute that h can be written as, 

-I(det J,r'^^^'Jl'p.P, - (det ^i)-V'^, (60) 

with h as in f l57|) . Concerning the function we want to match it with, we first of all need 
more info about A(„). Using the representation f|36|) of J, we can put 

Ji \ / 1 Jf Vi2 



from which it follows that 

n+l 

A(i)a'-^ = det( J + aPa) = (det Ji) det( Ja + a/„), (62) 

1=1 

and hence that A(„) = (det Ji)(tr J2). It follows from flS^ that 

(det Ji)"^//(„) = i g'^'pap, + i(tr J2) J^^p.p,- + (det Ji)-'W^n). (63) 

It looks by far not obvious that the expressions fl60p and (16 3 p would differ only by one term 
(up to irrelevant functions of time only). We shall compare them indirectly by computing 
their derivatives with respect to the x'^ and pa- The latter is easy and gives 

dh (9 

— = g''% + r%, — ((det J.r'H^n)) = 9"%. 

which is in line with the result we want to prove. For the other derivatives, the compu- 
tations can be written in a somewhat more compact form if we use the basis of vector 
fields 

adapted to the connection, rather than the coordinate derivatives with respect to x"". One 
can verify, recalling that = dtp'^/dx"', that Xa{h) can be written as 



+ ^-(det^O"V.^^. 

Making use of the scKt properties of J and Ji, plus the property dH) for Ji, this expression 
considerably simplifies and finally reduces to 

Xa{h) = \akg'^PiPa - laiiJl^'ikJaJlPiPj + \<^aJlPiPj 

+ (det JO yjl^. 
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where a as before stands for d{\x J\ Note in passing that, for example, J*-^ = ■ 
The computation of Xa ((det Ji)~^H(^n)) is much easier and gives 

((det Jl)-'H^n)) = \0laJl'P^pJ + (dct J^-^^^, 

with a = (i(tr J2). We need to make three more observations now. The first is that 

Xa{J^'pipb) = J\aPiPh = ^akg'^pipa, 

which takes account of the first term of Xa{h). Secondly, we recall the difference between 
aa and a^, as obtained in ( 15^ . which makes that the first term of Xa ((det Ji)~^if(„)) 
matches two terms of Xa{h). Finally, the terms not containing momenta also match as a 
result of (l42l) . taking into account that ^2 = —d2V. The conclusion now is that 

Xa{h) = Xa ((det Ji)-^i/(„) + J'%P,) . 

Remember that the two functions under consideration are regarded here as depending on 
the x", pa and time t (along solutions of the driving equations). Since their derivatives 
with respect to the x'^ and pa are the same, the conclusion is that they are indeed equal 
up to an additive function of time. □ 

The idea now is to try to get rid of the second term on the right in fl58|) by a further, suit- 
able canonical transformation. It is of some interest to look at this kind of question in all 
generality and to observe that it can be resolved indeed by a suitable point transformation, 
which necessarily must be time-dependent however. The notations used in discussing this 
general question below have nothing to do with any of the specific situations encountered 
so far. 

Lemma 3. Suppose that Hi{q,p, t) and H2{q,p, t) are two functions which differ by terms 
linear in thepi. Then, there exists a point transformation, {q,p) -H- {Q,P) say, such that 
the transformed Hamiltonian of the system with Hamiltonian Hi becomes the function H2 
expressed in the new variables. 

Proof. By assumption, we have Hi = H2+p^{q, t)pi for some functions p*. If F{q, P, t) = 
PiQ'^{q, t) is the generating function of an as yet unspecified point transformation, we know 
that the transformed Hamiltonian of the system with Hamiltonian Hi will be given by 



so that the desired effect requires that each of the satisfies 



dt dq^ 

In other words, we need n functionally independent first integrals of the equations = 
p^{q,t), which can be done in principle and is of course the same as saying that we have 
to integrate those equations. Note that even if the given would not depend on time, 
this procedure can only work with a time-dependent canonical transformation. □ 
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In the case of interest, we are looking with equation (!58|) at hnear terms of the form 



p%x,t)pa = r'{x,y{t))pi{t)p, 



We can try to find first integrals of the equations x"" = t) of the form u'^ = u'^{y{t),x), 
i.e. which are such that the time-dependence originates from solutions y{t) of the driving 
equations. This means that we have y'' = J'^^pi, and the first integral condition becomes 



We can now prove one of our main results, for which we go back to the supplementary 
assumption of the end of Section 4. The eigenfunctions u°'{y,x) which were introduced 
there will now be used as new coordinates for the driven system, along solutions y(t) of the 
driving system, and we will denote the corresponding conjugate momenta (for reference 
to the notations used in the Euclidean case in [TU]) by s^. 

Proposition 4. Assume that the equation det(J — XP2) = has n functionally inde- 
pendent solutions u°-{y,x). Then, the canonical transformation {x"',pa) O {x°',Pa) with 
generating function [54\ ), followed by the canonical transformation {x°',Pa) ^ [u"", Sa) with 
generating function F{x,s,t) = SaU°'{y{t),x), has the effect of transforming the Hamilto- 
nian h of the driven system into the function (det Ji)~^H(ny 

Proof. We know from Lemma 2 that the first step brings the Hamiltonian h into the 
form (!58|) . According to Lemma 3, the second step will eliminate the linear terms in 
the momenta in (!58l) . provided the functions u°'{y,x) have the property of making the 
right-hand side of (1M|) vanish. But the u°- satisfy the relations 0451) . from which it follows 
by a left action of Pi that 



At this point, it is important to be aware of another general and in fact very simple result. 

Lemma 4. Assume that a Hamiltonian system has a Hamiltonian K{q,p,t) of the form 
K = 'y(t)H{q,p,t) , whereby H is a first integral of the system and 7 is an arbitrary 
function oft only. Then H in fact cannot explicitly depend on time and the time- dependent 
Hamilton- Jacobi equation for K{q,p,t) reduces to the autonomous one for H . 

Proof. It is obvious that (with respect to the standard Poisson bracket) we have 
{K,H} = i{H,H} = 0, so that the property H = reduces to dH/dt = 0. The 
Hamilton- Jacobi equation for K then reads 




(64) 




(65) 



The desired result now follows by raising an index. 



□ 
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and looking for a complete solution S{q, t, ai) of the form S = W{q, ai) — «i / 'jdt imme- 
diately reduces it to the Hamilton- Jacobi equation 



Since det Ji, along solutions of the driving system, is a function of time only and Hi^n) 
is known to be a first integral of the driven system under the same circumstances, it is 
clear that the assumptions of Proposition 4 precisely bring us in a situation where we can 
draw the quite surprising conclusion that H(^n) will no longer be time-dependent in the 
{u, s) coordinates and separability of the Hamilton- Jacobi equation for the driven system 
is essentially a matter of separability of -ff(n). It remains to convince ourselves that the 
Hamilton- Jacobi equation for a system with as Hamiltonian is indeed separable. Now 
comes a rather subtle point in the argumentation. The point is this! We want to test 
separability of by using criteria which have an intrinsic, i.e. coordinate independent 
meaning. As such, it is the function which matters, expressed in any kind of canonical 
coordinates. Transformation formulas from one set of coordinates to another can depend 
on external parameters then, if needed, but should not be regarded as depending on time 
because time-dependent canonical transformations do more to the Hamiltonian function 
than just expressing it in the new variables. 

To be concrete now, the sufficient conditions for separability which we want to invoke are 
intrinsic indeed: they comprise the existence of a special conformal Killing tensor J for the 
kinetic energy metric of the Hamiltonian, plus a corresponding condition for admissible 
potentials V. The condition d{JdV) = in [2] for example, which actually corresponds to 
the particular case of the cofactor condition = —dW when the forces are conservative. 
It is well known (see e.g. [1]) that if the scKt involved in these conditions has functionally 
independent eigenf unctions, then these are separation coordinates for the Hamilton- Jacobi 
equation. It so happens that we have needed these ^^-coordinates already to prove in 
Proposition 4 that the original Hamiltonian h of the driven system can be transformed 
into the function H{n) (up to a factor). But having shown this way that Hamilton- 
Jacobi separability becomes a question about the function it is more appropriate to 
put this function to the test in the set of canonical coordinates (x",pa) which naturally 
presents itself prior to introducing the separation coordinates. The delicate issue alluded 
to above is that, when we subsequently want to pass to the new variables (m", Sa) again, 
the interpretation for this part of the story is that the functions u"'{y,x) are regarded 
then as depending parametrically on the y-coordinates of the driving system (not as a 
time-dependent transformation along solutions y{t) of that system). 

Going back to the expression f lB^ of -ff(n), taking into account that the function -f^{n+i) 
in dnn]) is actually a constant parameter (along solutions of the driving system) which we 
called in ( !T5|) . we have that 



This is the right expression for activating our test because all the ingredients we need for 
that have been prepared in Section 4. 




for the autonomous function H. 



□ 



H^r.) = |(det Ji) rtaPfc + (tr J2){E' - W\y)) + W^^){y, x). 



(66) 
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Proposition 5. J2 is a special conformal Killing tensor for the metric associated to the 
quadratic terms in ^6^) and the remaining terms satisfy the conditions for an admissible 
potential for Hamilton- Jacobi separability. 

Proof. We know from Proposition 1 tliat J2 is a scKt for g2 = {gab{x)) and recall that 
the characterizing property ([21) of scKts, when expressed in terms of the underlying type 
(1,1) tensor field is given by f|T2l) . with a = d{tTJ). It is then clear that the same J 
is also a scKt with respect to any constant multiple of g. In the case of the quadratic 
terms in f l66|l . we are precisely looking at a constant multiple of the metric g2 since 
(det Ji) is a function of the external y-parameters only, whence the conclusion about 
J2. We further observed above that the condition for admissible potentials is a reduced 
form of the cofactor condition that (cof J)fi should be closed. For the situation at hand 
we already know from Proposition 2 that the function Wi^n) satisfies this condition with 
respect to the scKt J2. It remains to show that the same is true for the remaining term 
in (166|) which in fact, since the factor {E^ — W^{y)) can be treated as a constant here, 
amounts to saying that the function tr J2 satisfies the condition. But this is trivially the 
case, because it follows from that for any tensor J with vanishing Nijenhuis torsion 
(cof J)d(tr J) = rf(det J). □ 

We sum up the main results about the driven system now. We know since Section 3 
that the driven system, along solutions of the driving system, has n integrals which 
are in involution with respect to two Poisson structures, one of which is the standard 
one when using coordinates adapted to the integrable distributions K and K-^. Under 
the assumption that the characteristic equation det (J — XP2) = has n functionally 
independent solutions u°'{y,x), we have seen in Proposition 4 that the given Hamiltonian 
h of the driven system, can be transformed into the function (det Ji)~^if(„). The fact 
that all the H(^i) are first integrals then implies that 

OH 

^ + (det Ji)-i{iJ(,),if(„)} = 0. 

But this, in view of the involutivity, actually means that none of the integrals will be 
time- dependent when expressed in the canonical coordinates {u"',Sa). This is in line with 
the results of Proposition 5 which mainly says that due to the existence of a scKt J2, 
the Hamilton- Jacobi equation of the (autonomous) function if(„) is separable. In fact, 
referring to the results which can be found in [1] for example, we can state more precisely 
that we are looking at a separable system of Stackel type and that the eigenfunctions 
u°'{y, x) of J2 are orthogonal separation coordinates. It should then be true indeed that in 
those coordinates we have n time-independent quadratic integrals in involution. Note that 
for the function H(^n) , for example, this time-independence means among other things (see 
the expression fl66l) ) that the function W(^n) — (tr J2)W^, when passing from the coordinates 
{y,x) to the coordinates {y,u{y,x)), should become independent of the ^/-variables. 

To give some further backing for all these rather subtle properties, we carry out more 
explicit calculations in the Appendix. The programme there is that we complete the 
recursive scheme started in Section 3 by proving, in analogy with the expressions (!22l 
[23l) . explicit results for the different parts of all the y4(j) tensors and computing explicit 
expressions for the corresponding functions 
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7 Examples 



For examples which exhibit all features of our generalization, we probably should think 
of systems with at least three degrees of freedom and a non-Euclidean metric. But going 
into the details of such applications would substantially add to the length of this already 
rather long paper, so we will keep this for a separate study and limit ourselves here to 
simpler two-dimensional situations. The first example is inspired by an integrability study 
in [T7] on what were called 'generalized Henon-Heiles systems', which are SODEs of the 
form 



Compared to older case studies of integrability of Henon-Heiles systems, the general- 
ization comes from the extra parameters a and m, which are motivated by allowing a 
Lagrangian description of the system in which the Hessian of the Lagrangian need not 
be (a constant multiple of) the unit matrix. The investigation carried out in [17] mainly 
consisted in looking for all possible parameter cases for which the system has two indepen- 
dent quadratic integrals. It led to the identification of three new cases, which are in some 
sense degenerate cases, because either a or m is zero, meaning that the Lagrangian one 
originally thought of is degenerate. As a result, there is no corresponding Hamiltonian 
which in the standard cases is always available as the first of two integrals in involution. 
A subsidiary question then was: to what extent can the two integrals in those degenerate 
cases be understood as being in involution, and it was argued that this can be resolved 
by constructing, in principle, a suitably adapted non-standard Poisson structure. As the 
equations in the case that either a or m is zero clearly exhibit partial decoupling, there is 
a good chance that those degenerate cases actually fit within our present theory, which is 
what we will discuss now. 

Consider the case where m = and & = 0, so that the system reduces to 



These equations are of course easy to solve without further ado, but they must serve here 
in the first place to illustrate various aspects of our theory. The second equation plays 
the role of driving equation. For consistency with the notations in the preceding sections 
we rename the variables as q2 = y, qi = x and write 



qi = -cigi + bql - aq^, 
h = -C2g2 - 2mgig2- 



-cigi - aqj, 



-C2g2- 



y 



-C2y, 

—C\X — ay^ . 



X 



Putting Py = y and p_ 



•x 



X the two quadratic integrals read 



^2 



1 



lr,2 _L V 7/2 

\{cx - 4:C2)pI + 2ayp^py - 2axpl + |ci(ci - 4c2)x^ 
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Obviously F2 is a Hamiltonian for the driving equation and is likely to be identifiable 
with the function H(^2) in the theory (see (1271) knowing that n = 1 here). The idea now 
is the following. Since Fi is a first integral of the complete system, its quadratic part 
identifies a Killing tensor A. Looking at the tensor J of which A is the cofactor, this 
may or may not be a scKt in general. It was shown in [H] however that this will always 
be the case for the Euclidean metric in dimension 2. If the component Ji is nonsingular 
therefore, we must be in a situation covered by our present theory and all the features we 
discussed should apply, with Fi = = H, the Hamiltonian of the quasi-Hamiltonian 
representation (l24l) . The extra assumption that the driven system should have a genuine 
potential, parametrically depending on the ?/- variables, is also automatically satisfied here 
by dimension. 

From the expression of Fi, we see that 

^=(-o'r '"iV (67) 



2ay c\ — 4c2 

and A = cof J, where J is the tensor with the following matrix components 



-2ay —Aax 

This tensor is indeed a scKt with respect to the Euclidean metric, and its Ji component 
is nonsingular (assuming ci 7^ 4c2), so we are in business and A^ = cof Ji = 1. For 
completeness, observe that 

/i = —C2ydy — (cix + ay'^)dx. 

in this example and that A^ = —dW indeed, with 

W = ici(ci - 4c2)x^ + a(ci - 2c2)xy^ + |aV- 

The function h for the standard Hamiltonian representation of the driven equation (with 
parameter y) is given by 

h = ^pI + ay'^x + |ciX^. 

Let us now further illustrate the subtleties of the theory, as explained in Sections 4 to 6. 
We start by computing the function h as defined in (l55l) . The tensor J2 here reads 

/ 4a\2 \ d 

J2 = - 4ax H — — ® dx. (69) 

\ ci — 4c2 / ox 



and the new momenta defined in (l52l 1531) become 

Py = (Cl - 4:C2)Py, 

Px = Px- 2aypy. 
In those new variables, the expression for Fi becomes 

Fi = i(ci-4c2)p^-2oV(ci-4c2)p;-2ax(ci-4c2)'pJ 



+ |ci(ci - 4c2)x^ + a(ci - 2c2)xy'^ + ha^y"^ 
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which is in agreement with (!63|) . For the computation of h on the other hand, we need the 
generating function (154|) of the time-dependent canonical transformation {x^Px) -H- 
We see from (168|) that Jl = —2ay = dip/dx with ip = —2axy. The generating function F 
thus reads 

F{x,px,t) = xpx - 2axy{t)py{t), 

and computing its partial time derivative involves making use of the driving equation. 
The resulting expression for h is found to be 

T 1 ~2 r\ 2 2 ~2 r\ / A \ ~2 r\ ~ ~ 2 12 *^^2 

h= 2Px + 2a y p - 2ax[ci - 4c2)p - 2aypxPy + ay x + ^c^x H — . 

Cl — 4C2 

One can verify that h and = Fx indeed verify the requirement (158|1 of Lemma 2 to 
within an additive function of time, which is the function 

4a y py 



2(ci - 4c2) 

and of course can be ignored in writing down Hamilton's equations. For the final canonical 
transformation to be applied to h, we need the eigenfunction u{y, x) of J2, which is found 
to be (see ([69])) 

4a^l/^ 

u{y, x) = 4ax. 

Cl - 4C2 

The time-dependent canonical transformation with generating function F{x, s,t) = s u{y(t),x) 
now transforms the relevant part of h into the function 



(det Ji)-iFi = 8a' + iw(ci - Ac2)pl + l^^^v' + ^2^'^ 

C\ — 4c2 oZa 



It then easily follows by taking into account that F2 is a constant along solutions of the 
driving equation, that 

Fi = 8(ci - 4c2)a^.^ + uF2 + '^^^2^''^ ^'^ 
which indeed no longer depends explicitly on time. 

The other two degenerate cases in [17] also fit within the present theory to some extent, 
but are even more peculiar. One can verify that the two integrals in those cases again can 
be understood by the fact that the system is a driven cofactor system. This also explains 
why the equations partially decouple. But it so happens that the Ji part of the scKt J 
in those cases is zero, so that we are not in the generic case of a nonsingular Ji, which 
was the assumption in the preceding sections. 

With a second simple example, we want to illustrate and test mainly the beginning of 
our theory, before suitable coordinates for partial decoupling have been identified, i.e. the 
situation covered by the conditions of Definition 2. Consider the system 

qi = 5gi - 4^2, 
q2 = qi- 
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The connection coefficients of the SODE connection are all zero and the (1, 1) tensor $ 
has the following matrix representation 



($}) 



df\ / -5 4 



dq^ J \-l 



With such a simple, constant Jacobi endomorphism, finding a distribution K which entails 
submersiveness of the system, i.e. which satisfies the conditions ([5]) of Definition 2, is 
simply a matter of looking for a 1-dimensional eigenspace of $. We choose 

( d d 
K = sp { — + 



dqi dq2 

Note that submersiveness now is ensured and is not related to the existence of a cofactor 
representation of our system. The problem of detecting such a representation is quite 
interesting in its own right. It bears to some extent resemblance to the inverse problem 
of the calculus of variations, because there is a certain freedom in selecting a multiplier 
matrix g ffist. For example, making the obvious choice of the unit matrix for g, or 
expressed in more mechanical terms, associating the given equations with a the standard 
kinetic energy T = + ^'2) "^i^ i^ot work! Indeed, one can verify that with this g, there 
is no special conformal Killing tensor J for which the right-hand sides of the equations 
will satisfy the condition Djfi = 0. On the other hand, the following choice for g turns 
out to be appropriate. Take 

then, a scKt J with respect to g is easily found to be 

( T ] = f '^'^^ ~ "^^^ ^1 + 
^ V 1i + 8q2 -4gi+40g2 

But more importantly, the nonconservative forces of our system are now of the appropriate 
form. That is to say, multiplying the right-hand sides of the equations with g, the 1-form 
/X is found to be 

fi = 4(gi - q2)dqi + (5gi + Aq2)dq2, 

and computing the cofactor tensor A of J, one can verify that indeed is closed (or 
equivalently -Dj/i = 0). We can now compute the orthogonal complement of the distri- 
bution K and see whether the final conditions for a driven cofactor system are verified. 
Since gu + (721 = 0, it easily follows that 



= sp 



_d_ 

dqi 



We have 

D"fi = Adqi ® dqi + 5dqi ® dq2 — 4:dq2 ® dqi + Adq2 ® dq2. 

It is clear then that fi{K^, K) ^ 0, whereas obviously d^{K,K) = by dimension. 
Hence, all requirements of Definition 2 are met. Integrating the distributions K and K-^, 
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suitable coordinates for the decoupling are found to be x = g2, y = Qi — (I2 and the 
transformed system becomes 



y = %, 
X = X + y. 

We leave it to the reader to verify from here, exactly as we did with the first example, 
that all the other features of our theory hold true, in particular those related to the 
consecutive canonical transformations relating the Hamiltonian h of the driven system to 
the quadratic integral -ff(i). 

8 Appendix 

Consider the recursive scheme for the tensors in particular the relations flTS]) for 
i = 1, . . . ,n. We have seen that taking i = n led to determining relations for three of 
the four parts of (see (l22l) ). while the remaining block A(^n)i had to be obtained from 
taking subsequently i = n — 1. This procedure can be continued all the way down, and 
the essential features of the recursion are captured in the following statement. 

Proposition 6. Suppose that we know the tensor parts /l(j+i)2i, A(i_|_i)i2 and y4(j+i)2 and 
that they satisfy the identities 

<^l^(j+l)12 + -^^12^(^+1)2 = = A(i_,.i)2l Jl + ^(i+l)2<^21, (70) 
J2lA(i+i)i2 + J2^(i+1)2 = A{i+l)2lJl2 + ^(j+l)2</2- (71) 

Then the following are determining equations for the completion of the construction of 
and for the next step in the recursion 

^(j+i)i = ^(i+i)-^r^ ~ 'Ji^Ji2A{i+i)2i, (72) 

A{i)21 = — <^21^(j+l)l — <^2^(j+l)21, (73) 
A{i)l2 = — ^(i+l)l<^12 — ^(i+l)12-^2, (74) 
A(i)2 = A(j+i)P2 — -^21^(i+l)12 — -^2^(i+l)2- (75) 

Moreover, the newly obtained parts of will satisfy the same three identities as those 
assumed for y4(i+i), while the completion of A(^i^i) will give rise to the following supple- 
mentary identity: 

JlA.(i+i)i + Jl2^(i+l)21 = + ^(i+l)12-^21- (76) 

Proof. Consider the recursion relation f|T8l) which in fact has two parts. The idea is to 
compose each of those parts on the left and on the right with one of the projectors Pj. 
This gives rise to a total of 8 equations. For example, acting with Pi on both sides of 
the relation JA(^i+i) + P2^(j) = A(t+i)-^Af (which we can call a (Pi, Pi) action for brevity) 
implies that we must have 

JlA(i+i)i + Jl2^(i+l)21 = A(j+i)Pi, 
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from which the determining equation (172|) follows. Likewise, a two-sided {P2, Pi) action 
will generate ( !73|) and the double P2 action generates (175|) . The remaining (Pi,P2) com- 
bination merely confirms the first of the assumed identities ( !70|) . Starting from the other 
part in (1181) . the determining equation (I74p will follow from a (Pi, P2) action. The (P2, Pi) 
combination confirms the second of the assumed identities in (l7Up . The double P2 action 
gives rise to another determining equation for A(j)2, which is consistent with the first one in 
view of the identity ( I7T]) . Finally the double Pi action implies that, again for consistency, 
y4(j+i) must further satisfy the identity fl76|l . It is then a straightforward computation to 
verify that the obtained blocks A(j)2i, A{j)i2 and A(i)2 will satisfy corresponding identities 
of the form ( 170|) and ( 1711) in view of those assumed for A(j+i). This concludes the full 
recursion step. □ 

The identities which the different parts of all tensors satisfy are important to get to 
a considerable simplification of the recursive scheme. Indeed, as will be shown now, it 
turns out that knowledge of the block ^(i)2 suffices to determine the three other blocks of 
Moreover we can set up a recursive scheme to determine ^(j)2 from A(j_|_i)2 and we 
shall see that this procedure brings the tensor J2 back into the spotlights. 

It follows from the identities fiTOl) that for each i = 1, . . . ,n: 



A[i)12 = —Jl ^Jl2^(i)2, (77) 
^(i)21 = — ^(i)2-^21<>^r\ (78) 



and subsequently from the defining relation ( 1721) that 

A(,)i = A(i) Jf ^ + JrVi2A(i)2J21 Jf'. (79) 

Making use of ( 1771 1781) in the identity (17T1) it is easy to see that this expresses the com- 
mutativity 

J2v4(i)2 = v4(,)2J2. (80) 

Finally, the determining equation ( 175|) reduces to 

^(j)2 = ^{i.+l)P2 — </2^(i+l)2- (81) 

Lemma 5. The blocks ^(i)2 = P2 ° ^(i) ° P2 of the A(^i-^ tensors which determine the 
quadratic part of the first integrals H{i) are recursively given by 

n—i 

A(,)2 = A(,+i)P2 + 5^(-l)^A(,+,+i) J2^ i = l,...,n-l, (82) 
and all other parts 0/ follow from A(j)2- 

Proof. We know that A(„)2 = A(„+i)P2, with A(„+i) = det Ji. The recursive relation 
(!82|) then easily follows from ( IHT]) by induction. The last part of the statement has already 
been proved above. □ 

Obviously, we are now in a position to venture computing a more explicit expression for 
the functions and it turns out that it is most appropriate to do this in terms of the 
momenta p again. 
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Proposition 7. The quadratic integrals H(^i) (for i = 1, . . . ,n) are given by 

H^,) = ^Al^PaPb + |A(,) Jf M« + (83) 

where 

n—i 

= A(,+i)^'^'' + 5^(-l)^A(,+,+i)(J2^r'- (84) 
i=i 

Proof. Recall that H{n) has already been computed (see (156!) ) and is indeed of the form 
(!83|) . From (126|) . we further have that 

Observe first that raising indices in ( IHTl) gives rise to the formula 

Aab A 1 b Ada 

and it subsequently follows from ( 1771) and ( 1791) that 

We now compute each of the quadratic parts of in terms of the p, using fl52 l 153|) . The 
first term, for example, becomes 

It SO happens that only the first term on the right in this expression survives, in other 
words all the other terms in the end cancel out when we proceed in the same way with 
the two other quadratic terms of -f^(j)- The formula f l83|) then readily follows, while f lH^ 
merely is the contravariant form of fl82|) . and will turn out to be useful further on. □ 

As was the case with -ff(n), we recognize in the expression for the other H(j_) part of the 
constant E"^ of the driving system. Explicitly, with the help of ( lT5l) . it is clear that along 
solutions of the driving system, the -ff(j) can be written as 

= iAgp.p, + -^^{E' - + ly^, (85) 

whereby we recall that det Ji = A(„_|_i). We shall now finally illustrate that applying the 
parameter-dependent coordinate change (x'^,Pa) O (M"(x,y), Sa) turns the expressions for 
the into functions which no longer depend on the parameters ?/*, and hence, in their 
interpretation of first integrals of the driven system along solutions of the driving system, 
become effectively time-independent quadratic integrals. 

According to Proposition 3, the u"" are eigenfunctions of J21 so that the coordinate ex- 
pression of J2, in the variables {y.u) takes the simple form 

d 

OU"' 
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Going back to the result (162|) . it follows that 



n+l n 

A(i)a^-^ = (det Ji) det( J2 + a/„) = (det Ji) + a), 

1=1 b=l 

which in turn implies that 

A(i+i) = (det Ji)a„_i(u), i = 0,...,n, (86) 
where (Tj{u) denotes the elementary symmetric functions. 

For information about the functional dependence of the other terms in the potential part 
of the we go back to the double cofactor representation of the overall system, which 
was the start of the recursive scheme in Section 3. We know that with A{a) representing 
the cofactor tensor of J + aP2, the force terms /i of the overall system satisfy the relation 
A{a)fi = —dW{a), or equivalently 

det( J + aP2)/i = -{J + aP2)dW{a). 

Projecting this relation under Pi it follows that 

det(J + aP2)/ii = -JiPi{dW{a)) - J2iP2{dW (a)) . (87) 

It is straightforward to compute that under the transformation {y,x) -H- {y,u{x,y)) the 
tensors involved acquire the following coordinate expression 

(d du^ d \ ■ d du"" 

■7^ + ^^] ®dy'^ P2 = ^^du" - ——®dy^ , 



and 



Ji = J]{— + — — 1 ® dy\ J21 = J-—— ® dy' 



^dy' dy^ du^ J ' '' dx"- du^ 

Concerning the left-hand side of ( IHTl) we have to remember also that (det Ji)/ii = —dW^, 
since the driving system is of cofactor type with scKt Ji. It then readily follows, using 
the property fl65l) . that the coordinate expression of fl87|l reduces to 

det{J2 + aP2) ^ ' 



Qyk Qyk ' 

and this for all a. The expansion (162|) thus implies that the functions — (A(j)/ det Ji)W^ 
do not depend on the y parameters, which together with (186|) confirms our objective for 
the potential part in H(iy 

It remains to look at the terms quadratic in the momenta. Expressed in the momenta 
associated to the w-variables, they become 
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We know that in particular A^^^ = (det Ji)g°'^ (cf. (156|1 ). Since J2 is a scKt with respect 
to g2-, it is symmetric in its covariant or contravariant representation. But J2 is diagonal 
in the u coordinates, hence the same is true for the transformed g2- Let us rely here 
further on the fact that we proved by indirect means, mainly as a result of the simple 
Lemma 4, that will be time-independent when expressed in the (u, s) variables. 
The net conclusion then is that the diagonal elements of the transformed g2 must be the 
product of (det Ji)~^(?/) with a function depending on the u variables only. It subsequently 
follows from the explicit expression of A'^^-^ in f lH^ and the fact that the ratios A(j)/ det Ji 
also are functions of the m- variables only (see ( l86l) ). that the quadratic part of all 
will indeed become time-independent as well, when these functions are looked at as first 
integrals of the driven system along solutions of the driving system. 

To complete the picture it is perhaps worth repeating that the final in that hierarchy 
is in fact the Hamiltonian H in the quasi-Hamiltonian representation of the full system 
and that (det Ji)~^A(i)2 is the cofactor of J2. Acting with J2 on the expression ( l82l) for 
i = 1, we thus obtain a polynomial expression satisfied by J2- It is then comforting for the 
internal consistency of our results that one can verify that this identity indeed expresses 
the Cayley-Hamilton theorem applied to J2. 
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